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ABSTRACT 

We perform TV-body simulations of young triple systems consisting of two low-mass 
companions orbiting around a significantly more massive primary. We find that, when 
the orbits of the companions are coplanar and not too widely separated, the chance of a 
collision between the two companions can be as high as 20 per cent. Collisions between 
one of the companions (always the less massive) and the primary can also occur in 
systems with unequal-mass companions. The chance of collisions is a few per cent in 
systems with more realistic initial conditions, such as with slightly non-coplanar orbits 
and unequal-mass companions. If the companions start widely separated then collision 
are very rare except in some cases when the total mass of the companions is large. We 
suggest that collisions between members of young multiple systems may explain some 
unusual young multiple systems such as apparently non-coeval companions. 

Key words: methods: iV-body simulations - binaries: general - stars: formation. 



1 INTRODUCTION 

It appears that m any, if not mo st , stars form in mul- 
tiple systems (see iMathieul 1 1994 iGoodwin et al.l 120071 : 
lGoodwinll2010h . Indeed, an increasing proportion of both 
pre-main sequence (PMS) and main sequence (MS) mul- 
tiple systems are being found in triple or higher-order 
systems (e.g. ICorreia et al.l 120061; iTokovinin et alj l200rj ; 



lEggleton fc TokovininlbOOSHLaw et al.ll2010j) . There is also 
evidence for dynamical decay in young stellar systems sug- 
gesting that the initial triple and hi gher-order fraction 
may be even higher than the field IHaisch et alj 120041 : 



IConnellev. Reipurth. fc Tokunagall2008al lbl). 

Observations have found a number of young stars with 
unusual companions, that is companions w i th unexpected 
colours llHartigan. Strom, fc Stroml 1 1994 ; Koresko et al.l 



19971 : iDuchene. Ghez. fc McCadell200"3t Hartigan fc Kenvori 
20031 : iPrato. Greene, fc Simor] 120031 : |k raus fc Hillenbrand! 
2009). The interpretation of these unusual compan- 



ions is complex with some authors claiming some non- 
coevality between components in multiple systems (e.g . 



Hartigan. Strom, fc Stroml 19941: lHartigan fc Kenvonl 12003 
Prato. Greene, fc Simonl 120031 ; iKraus fc Hillenbrand! 12009^ 



although at least some such objects may be the result of an 
enshroude d companion as appears to be the case in T Tauri 
itself (see iDvck. Simon, fc Zuckermanl Il982 ; Ratzk a~et al.l 
2009). Non-coevality is suggested by the two objects having 
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positions on PMS track s that cannot be ex plained by dif- 
ferent masses. However, iBaraffe et al. | l|2009h has suggested 
that apparent non-coevality may be due to differences in the 
accretion history of the protostars leading to different evo- 
lution along PMS tracks. Within a multiple system this may 
be due to different accretion histories (through the different 
angular momentum of infalling material?), but we suggest 
that it is possible that it could be due to a collision be- 
tween low-mass companions early in the history of the sys- 
tem. Indeed, a collision has been suggested to explain the 
underluminosity of the lo w-mass brown dwarf 2M1207B by 
iMamaiek fc Meverl <|2007h but described as 'improbable', we 
will show that such a solution might not be as improbable 
as one might first think. 

We perform iV-body simulations of a large ensemble of 
triple systems to investigate the frequencies of collisions in 
such systems. We describe the details of our simulations in 
Sec. [2] The results are presented in Sec. [3] and followed by 
discussion in Sec. [4] 



2 SIMULATIONS 
2.1 Initial conditions 

We set up a three-body system of young stars in which two 
low-mass companions of masses M2 and M3 are orbiting 
around a primary star of (a higher) mass M\ . The compan- 
ions are assumed to form with in a circumstellar disc (e.g. 
Sta matellos fc Whitworthll2u09h . We simulate the evolution 
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of the system once the majority of the disc has disappeared 
(by accretion or removed by feedback?) and so we can de- 
scribe the evolution as a simple N = 3 JV-body problem. It is 
quite possible that the companions interact or collide with 
each other whilst a massive disc is present, and they may 
well interact with the disc in such a way that one or both 
migrate inwards or outwards. Such situations are extremely 
interesting but require modelling the hydrodynamics of the 
disc. 

The system can be characterised by three parameters: 



2.1.1 Protostellar masses 

We simulate typical T Tauri systems with primary masses 
Mi = 1 or 2Mq . The total masses of the companions, M2 
and M3, have a range M2 + Ms = 0.1 — 0.6M@ in steps of 
O.O5M , and mass ratios of M2/M3 = 0.25, 0.5 and 1. This 
gives us 33 combinations of companion masses for each pri- 
mary mass. 



2. 1.2 Protostellar radii 

The radius of a young star in the systems is estimated 
and scaled from the mass-radius relation of low-mass main- 
sequence stars. In t his work, we use the empiri cal mass- 
radius relation from ICaillault fc Patterson! i|l990l ) which is 
for stars of mass ~ 0.1 — 0.5Mq in the Solar neighbourhood. 
We assume that the relation extends to very low-mass stars 
of a few tens of Jupiter mass. Since the PMS stars are young, 
their radius could be larger by a factor of a > 1. The radius 
R of a PMS star of mass M can then be written as 



R = 0.92a 



M 

M© 



R 



©• 



(1) 



The values of a used in the simulations are selected from 1 
up to 20: in steps of 0.5 for a = 1 — 10, and in steps of 1 for 
0=10-20. 



2.1.3 Initial positions 

Observations of PMS binaries show that they have a range 
of separations from sub-AU to a few hundreds of AU (e.g . 
Mathieu 1994; [Patience et al.l [20021 : iKonopackv et al. I l2007t 
Goodwinl 20101 ) . The peak in separations of PMS binaries ap- 
pears to be a round 100 — 200 AU with an excess above th e 
field (see e.g. IPatience et"aill2002l : IKonopackv et alj|2007l ). 
To cover the bulk of the separation range of PMS com- 
panions we perform simulations with companions between 
20 - 100 AU, 100 - 200 AU, and 200 - 300 AU. 

The two companions are placed randomly in the given 
ranges measured from the primary and usually both in the 
same plane (we also perform simulations to test the effect of 
non-coplanar motion by giving the third companion a small 
velocity component in ^-direction). The angular separation 
between the companions is forced to be greater than 45° in 
order to avoid initial states which result in a collision or ejec- 
tion almost immediately. They are given the correct velocity 
for a circular orbit, and both orbit in the same direction. 

For each value of Mi, Ms, M3, and a we perform an 
ensemble of 5 x 10 3 or 10 4 simulations. 



2.2 Numerical method 

The integrator that we use is a variable stepsize forth- 
order Adams-Bashforth-Moulton predictor-corre ctor (e.g. 
iBinnev fc Tremaind I2OO8I : iMathews fc Finkll2004h . This in- 
tegrator requires only two force evaluations per timestep to 
provide a solution with high accuracy. The fractional energy 
error of the calculations is kept below 10 -5 , and a system 
with unacceptable error is reintegrated from the beginning 
with a higher accuracy. If an error persists the result is omit- 
ted from analyses (reducing the number of the systems in 
an ensemble). 



2.3 Termination criteria 

Simulations can result in either a collision, an ejection, or a 
stable system. 



2.3.1 Ejections 

We consider a companion 'ejected' if it travels further than 
5000 AU from the primary. Whilst this might not always 
be unbound, we consider that a companion on an orbit of 
at least 5000 AU apastron distance will be very soft and 
likely to be perturbed by interactions with other stars in 
the cluster in which we assume our young multiple system 
has formed. 



2.3.2 Collisions 

Stars collide if the mutual separation is less than the sum 
of their radii as calculated above. There are two possible 
types of collisions that can occur: (i) collisions between com- 
panions (companion-companion collisions, or CCCs), and 
(ii) collisions between a companion and the primary star 
(companion-primary collisions, or CPCs). 



2.3.3 Stable systems 

We consider a system stable if there are no ejections or col- 
lisions within 5 Myr. This length of time covers the typical 
ages of clusters containing T Tauri systems. (As we shall 
see, stable systems are rare and so the exact length of this 
criterion is unimportant as long as it is > 1 Myr) . 



3 RESULTS 

Our fiducial system has a primary mass Mi = 1M© , and 
companion separations between 100 — 200 AU about the 
peak of the PMS separation distribution. The main result 
of the fiducial system simulations is as would be expected: 
most of the systems decay by ejection within ~ 100 cross- 
ing t imes, preferent i ally ejecting their least-m assive member 
(e.g. lAnosovalH986l : ISterzik fc Durisenlll998h . However, we 
also often find a significant number of companion-companion 
collisions (CCCs) but very few stable systems or companion- 
primary collisions (CPCs) . In ensembles with different initial 
conditions stable systems and CPCs can become important, 
however dynamical decay is most often by far the dominant 
outcome (see below). 
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Figure 1. The effect of changing the radius factor (a) in our 
fiducial system (with primary mass IMq and orbits between 
100 — 200 AU) on the frequencies of ejections (blue lines with di- 
amond symbols), CCCs (red lines with circle symbols), CPCs 
(purple lines with triangle symbols), and stable systems (cyan 
lines with square symbols). The frequencies are plotted against 
total companion mass M2 + M3 of the ensembles. In all cases 
M2 = M3. Increasing line thicknesses correspond to increasing ra- 
dius factors a = 1, 5, 10, 15, 20. Note that the fractions of CPCs 
and stable systems are negligible. 



In a three-body system such as the ones we are simulat- 
ing there is one major body (the primary) which dominates 
the potential. At sufficiently low-mass companions would 
barely feel each-other's gravity and would evolve as a plane- 
tary system (Mi 3> M2 + M3). As the mass of the compan- 
ions increases they will perturb each-other's orbits, but al- 
ways in the global potential of the primary (Mi > M2 + M3). 
As their orbits are perturbed by each other the companions 
will undergo a number of close encounters which will trans- 
fer energy from one body to the other. Often an encounter 
will provide one of the companions with enough energy to 
escape the system altogether - an ejection. However, a not 
insignificant fraction of encounters will be close enough to 
cause a collision. This collision is almost always between the 
two minor bodies in the system as they are the bodies whose 
orbits are most perturbed by the other. 



3.1 Companion-companion collisions 

In Fig. [1] we show the change in the number of ejections 
(lines with diamond symbols) and collisions (lines with cir- 
cle symbols) for our fiducial system with increasing total 
companion mass where both companions are of the same 
mass (M2 = M3). The increasing thickness of the lines shows 
the effect of increasing the radii of all the stars by a factor 
a = l,5,10,15,20 (moving upwards with increasing thick- 
ness to higher a for CCC, and downwards for ejections). 

Unsurprisingly, larger stellar radii result in more col- 
lisions, from a few per cent when a — 1 (the MS radius), 
to around 20 per cent when a = 20 (a very large PMS star). 
Collision timescales are similar to those of ejections and usu- 
ally occur within ~ 100 crossing times. 

The fraction of CCCs stays relatively constant for any 
given a as M2 + M3 changes from 0.1 to O.6M0 (from 10 




Figure 2. The probability of an encounter at a companion sep- 
aration ^r23, P{r23), as measured by the frequency of CCCs 
with increasing a. In all cases M\ = IMq and M2 + M3 = 0.3Mq . 
Black circles are for systems with mass ratio M2/M3 = 1 and red 
circles for M2/M3 = 0.25. Solid lines are the fits from Eqn. (J2)l for 
all data, whilst dashed lines show the fits for only a = 1 — 10. 



to 60 per cent of the primary mass). This is because even 
at 0.6Mq, the companions are still not the dominant con- 
tributor to the potential. And even though more massive 
companions require a greater escape energy (as they are 
more massive and escaping from a more massive system), 
the energy change in an encounter between more massive 
companions is greater by the same order. 

The number of CCCs changes in an unusual way with 
increasing a. The number of CCCs increases quite rapidly 
as a is increased from 1 to 5 to 10, but the relative increase 
in the number of collisions as a is increased from 10 to 15 
to 20 is small. 

A CCC occurs when the separation of the companions 
r23 is equal to the sum of their radii R2 + R3 ■ Therefore, 
the number of CCCs is the probability that the separa- 
tions become less than r23, -P(»"23). Fig. shows the in- 
crease in P(r23) with r23 for two example systems with 
M 2 + M 3 = 0.3M Q and mass ratios M2/M3 =0.25 and 1. 

The trends in Fig. [5] can be fitted very well with the 
Levy cumulative distribution function (CDF) of the form 



P (J23) = aerfc 



(V & A'23 



+ C, 



(2) 



where a, b and c are constants, and erfc is the complemen- 
tary error function. This CDF provides the best fit with 
meaningful properties, namely, dP(r) /dr = at r = and 
lim P(r) = constant. The values of the constants and co- 

r — ¥ 00 

efficient of determination (1Z 2 ) corresponding to each data 
set are summarized in Table[T] The solid lines in Fig.[2]show 
the best fits to high-a, while the dashed lines show the best 
fits to low-a (see below). 

That the data are well fitted by the Levy distribution 
can be explained by assuming the following, (i) As the dis- 
tance between the companions at collision is far smaller than 
the distance to the primary the encounter between the two 
companions is approximately a two-body problem, (ii) The 
effects of a rotating frame of reference about the primary 
are small at low r23. The problem can then be simplified to 
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Table 1. The constants in Eqn. J2]| obtained from nonlinear re- 
gressions of the data in Fig. [2] The goodness of fit is represented 
by the coefficient of determination TZ 2 (closer to 1 is better) . 



a 


M2/M3 


a 


b 


c 


V? 


1-20 


1.00 
0.25 


0.254 
0.207 


1.242 
1.265 


0.091 
0.079 


0.9971 
0.9945 


1-10 


1.00 
0.25 


0.232 
0.180 


0.972 
0.847 


0.084 
0.070 


0.9988 
0.9984 



one body of a reduced mass it, = M2M3 / (M2 + M3) orbiting 
around a fixed central mass of M = M2 + M3 . Let us now 
consider T23 as the separation between the masses fj, and M 
at the pericentre of the orbit. At this turning point, as the ra- 
dial velocity is zero, we assume further that (iii) the tangen- 
tial velocity (vt) of the mass li follows the one-dimensional 
Maxwell-Boltzmann distribution, i.e. f(v t ) oc exp(— v 2 /2a 2 ), 
where a is the velocity dispersion. The probability of the 
mass fx having vt from that for a circular orbit to hyperbolic 
orbits around the fixed mass M may be written as 



P(vt ^ Vciv) oc 



/2a 



dv t oc erfc(ti c i r /\/2<T) . 



Substituting v c - lr — y GM/r23 in Eqn. ©, we have 
P{r23) = aerfcfy / &7r23j , 



(3) 



(4) 



where a and b are constants. The constant c in Eqn. ((2} is 
due to 'head-on' collisions in which v t is ~ 0. 

In addition to the solid lines in Fig. [2] the dashed lines 
show the fits of the data with a ^ 10. Although these curves 
are formally better fits (higher TZ 2 , as shown in Table [TJ 
they tend to diverge from the data at higher a. These diver- 
gences probably indicate some complicated dynamics that 
we have not included in our derivation above (i.e. the effects 
of motion in a rotating frame of reference). 

The differential of P(r23) with respect to T23, p(r2s), 
is the probability density function (PDF) of the encounter 
having a separation at pericentre of T23, 



P(f23) 



b e- b/r23 
3/2 ■ 



(5) 



3.2 Systems of unequal-mass companions 

The companion mass ratio (M2/M3) plays an important role 
in the energy redistribution during close encounters between 
the companions. In iV-body systems, the objects in the sys- 
tems tend to distribute kinetic energy equally (equiparti- 
tion) . The systems thus usually consist of slow-moving-high- 
mass and fast-moving-low-mass objects. 

In our three-body system, most of the time only 
the companions are closely interacting with each other 
(see above). For systems with unequal- mass companions 
(M2/M3 <l), the equipartition of the kinetic energy usu- 
ally causes the lower-mass companion to be ejected from the 
system. For systems with equal-mass companions the com- 
panions need more close encounters before energy exchange 
is able to eject one object. The chance of collisions in the 
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Figure 3. The effect of changing the companion mass ratio 
(M2/M3) on the frequencies of CCCs (red lines with circle sym- 
bols), CPCs (purple lines with triangle symbols), and stable sys- 
tems (cyan lines with square symbols). Different line styles rep- 
resent different mass ratios: M2/Af3=0.25 (dotted lines), 0.5 
(dashed lines), and 1.0 (solid lines). In all cases Mi = IMq and 
a = 5. 



equal-mass systems is therefore higher than in unequal-mass 
systems. In Fig. [3] we can see that the number of collisions 
between companions (marked with circles) increases with 
the increasing mass ratio from M2/M3 = 0.25 (the dotted 
line) to M 2 /M 3 = 1 (the solid line). 

Figure. [3] also shows an interesting feature of sys- 
tems with unequal-mass companions: increasing numbers of 
companion-primary collisions (CPCs). We find that it is al- 
most always the lower-mass companion that collides with 
the primary star. This occurs after the lower-mass com- 
panion has been regularly perturbed by the more-massive 
companion orbiting the primary with a larger orbit. Some 
angular momentum is extracted from the lower-mass com- 
panion during each close encounter. This causes the orbit of 
the lower-mass companion to become more and more eccen- 
tric until eventually colliding with the primary star. CPCs 
are almost never seen with equal-mass companions. 



3.3 The effect of other parameters 

Apart from the radius factor (a) and the companion mass 
ratios, three other parameters can also affect the results. 
Firstly, the primary mass and the separation from the pri- 
mary to the companions (essentially the potential energy of 
the system). Secondly, the distance between the companions. 
And thirdly, the coplanarity of the orbits. 



3.3.1 Potential energy 

The potential energy of the system is set by the primary 
mass and the distance of the companions from the primary. 
The deeper the potential well, the more difficult it is for an 
ejection to occur (the more energy needed to be transferred). 
Therefore, the deeper the potential well the more encoun- 
ters are needed before an ejection occurs, and the greater the 
chance of a collision or the 'partial' ejection of a companion 
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Figure 4. The effect of changing the companion initial orbital 
range on the frequencies of CCCs (red lines with circle symbols), 
CPCs (purple lines with triangle symbols), and stable systems 
(cyan lines with square symbols). Different line styles represent 
different ranges of the orbital radii: 20 — 100 AU (dotted lines), 
100 - 200 AU (dashed lines), and 200 - 300 AU (solid lines). In 
all cases Mi = IMq , M 2 /M 3 = 1, and a = 5. 
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Figure 5. The effect of large initial separations between the com- 
panions (r2i) on the frequencies of ejections (blue lines with di- 
amond symbols), CCCs (red lines with circle symbols), CPCs 
(purple lines with triangle symbols), and stable systems (cyan 
lines with square symbols). Dotted lines represent the systems 
with M2/M3 =0.5 (the less massive in the smaller orbit), dashed 
lines with M2/M3 = 1, and solid lines with M2/M3 =2 (the more 
massive in the smaller orbit). In all cases Mi = IMq and a = 5. 



to a wider stable orbit. In Fig. [4] we decrease the separa- 
tion range of companions from 200 — 300 AU (solid lines), to 
100 - 200 AU (dashed lines), to 20 - 100 AU (dotted lines). 
The number of CCCs (red lines with circles) increases with 
decreasing separation (increasing potential energy). In the 
case of 20 — 100 AU separations the number of stable sys- 
tems is also significant (dotted cyan line with squares), in- 
deed, at low companion masses stable systems outnumber 
CCCs. Similar effects are also seen with increasing primary 
mass (not illustrated). 

3.3.2 Separation between the companions 

Initial separation is an important parameter. Frequent in- 
teractions leading to ejections and collisions occur far less 
often in systems in which the companions start widely 
separated. We demonstrate this with three ensembles 
whose companion's initial orbits are each restricted to be 
within two widely-separated ranges of r 2 = 50 — 60 AU and 
r 3 = 200 - 210 AU: (i) M2/M3 = 0.5 (the less massive in the 
closer orbit), (ii) M 2 /M 3 = 1, and (iii) M 2 /M 3 = 2 (the more 
massive in the closer orbit). 

We can see from Fig. [5] that most systems that start 
widely separated are stable. Systems tend to be more sta- 
ble if the lower-mass companion has the closer orbit, or if 
M 2 /M 3 <1; compare the dotted cyan line (M 2 /M 3 = 0.5) 
with the dashed cyan line (M2/M3 — 1) and the solid cyan 
line (M2 /M3 — 2) in Fig. [5] Similarly to our fiducial systems, 
CPC apparently occurs only in the systems with M2 /M 3 < 1 
(the dotted purple line with triangle symbols). 

Although the number of CCCs drops significantly, they 
do not disappear entirely except for very low companion 
masses. It is unclear what the initial architecture of young 
triple systems would most often be. Disc fragmentation 
would be expected to occur a t large radii (> 100 AU, 
IStamatellos fc Whitworthll2009h . however migration within 
the disc may rapidly separate companions. 



3.3.3 Coplanarity of the orbits 

To test the effect of non-coplanarity we perform some en- 
sembles of simulations in which the companion M3 is given 
a small velocity component in the z-direction just enough 
to make its orbit initially inclined by around (a) ~ 5° and 
(b) ~ 10° . We find that the number of collisions drops sig- 
nificantly, as shown in Fig. [5] from ~ 15 per cent at zero 
inclination (solid line at top), to a <5 per cent at 5° and 
10° inclinations (dashed and dotted lines). Non-coplanar in- 
teractions produce fewer collisions due to the introduction 
of a third dimension, however there are still a non-negligible 
number of collisions. 



4 DISCUSSION 

We have studied the dynamical evolution of coplanar triple 
systems of young stars. We find that collisions between mem- 
bers are not very unusual and can occur up to 20 per cent 
of the time in coplanar systems with two equal-mass com- 
panions. Collisions are mostly between the two companions 
rather than between a companion and the primary. We un- 
surprisingly find that stars with larger radii are more likely 
to collide, but that collisions are not uncommon at stellar 
radii only a few times the main sequence radius. 

We find that in more realistic situations with different 
mass companions and slightly non-coplanar orbits the num- 
ber of collisions is significantly lower, but can still be a few 
per cent. 

We conclude that collisions in young triple systems 
would not usually occur, but that they might hap- 
pen often enough that they could explain some un- 
usual systems such as some of the apparently non- 
coeval T Ta u ri sy st ems that have been observed 
l|Koresko et al.l Il997i ; iDuchene. Chez, fc McCadel 120031 : 
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Figure 6. The effect of changing the orbital coplanarity on the 
frequencies of CCCs (red lines with circle symbols), CPCs (purple 
lines with triangle symbols), and stable systems (cyan lines with 
square symbols). Different line styles represent different degrees 
of orbital inclination between the two companions: coplanar or- 
bits (solid lines), ~ 5° inclined orbits (dashed lines), and ~ 10° 
inclined orbits (dotted lines). In all cases M\ = IMq, M2/M3 = 1, 
and a = 5. 
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In our simple JV-body simulations we are unable to de- 
termine what would happen after a collision and what the 
product or products might be. Many collisions may be glanc- 
ing collisions stripping material off of one or both objects, 
or collisions m ay cause the objects t o merge to make a new, 
larger, object i|Lavcock fc Sillsll2005l ). 

Whatever happens, a collision would be expected to 
change the mass of the final object (s) either stripping ma- 
terial or forming a new object. Such rapid and violent mass 
changes would be expected to change the structure of the 
objects and cause them to be in unexpected p laces for their 
ages on the HR diagram ijBaraffe et al.l 120091 ) . Stripped or 
ejected material could persist around one or other object 
causing it to have a different (and unusual?) extinction to 
other ob jects in the system (an explanation for T Tauri it- 
self? See lRatzka et al'j|200gh . 

Without detailed hydrodynamic simulations and ra- 
diative transfer modelling it is impossible to know what 
the colours of the new obje ct (s) after collision might be. 
iKraus fc Hillenbrand! (|2009h find non-coevality with the 
most massive star appearing younger. If the mass determi- 
nations are accurate then this might suggest CPCs above 
CCCs (despite CCCs being more common in most of our 
simulations). However, the collision product will presum- 
ably be larger and so appear more massive, but it may 
also be hotter and so appear older? The collision product 
will presumably follow a different path towards the main se- 
quence than it would previously have done, not appearing 
the 'correct' colour and magnitude until it reaches the main 
sequence. Without detailed modelling it is very difficult to 
guess how a PMS collision product would look or evolve, but 
it seems reasonable to assume that it would very different 
to a 'normal' PMS star of the same (post-collision) mass. 

In summary, collisions would be expected in some young 



triple systems. They are not common, but are not so rare as 
to b e ignored in attempting to explain some unusual systems 
(e.g. iMamaiek fc Meverj|2007l ). 
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